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(N . Abstract 
> 

The predictions for the shell structure of metal clusters of the three-dimensional q- 
deformed harmonic oscillator (3D g-HO), utilizing techniques of quantum groups and hav- 
ing the symmetry u g (3)Dso 9 (3), are compared to the restrictions imposed by the periodic 
orbit theory of Balian and Bloch, of electrons moving in a spherical cavity. It is shown that 
agreement between the predictions of the two models is established through the introduc- 
tion of an additional term to the Hamiltonian of the 3D g-HO, which does not influence 
the predictions for supershells. This term preserves the u 9 (3)Dso g (3) symmetry, while in 
addition it can be derived through a variational procedure, analogous to the one leading 
from the usual harmonic oscillator to the Morse oscillator by introducing the concept of 
the Variable Frequency Oscillator (VFO). 
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C^. 1 Introduction 
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Algebraic models are popular in several branches of physics In ideal cases, algebraic 

models are able to describe in good approximation several properties of a physical system 
using a limited number of appropriate (usually collective) degrees of freedom. In addition 
to providing successful predictions by themselves, exactly soluble algebraic models can also 
serve as a useful testground for more sophisticated microscopic theories, in which heavy 
numerical work is inevitable. 

An algebraic approach to the structure of metal clusters E] has been introduced 
recently [SI El HI El Ej, taking advantage of the three-dimensional g-deformed harmonic 
oscillator (3D g-HO') |R)jL w hich is constructed using the techniques of quantum algebras 
(quantum groups) [se e flT| a nd references therein]. The symmetry of this oscillator is 
Uq(3)DsOq(3) CSl iM p lH l HH], while its derivation involves use of irreducible tensor 
operators under so 9 (3) E3 CHI- I n terms of the 3D g-HO a good description of the magic 
numbers of alkali clusters (up to 1500 atoms, which is the limit of validity for theories based 
on the filling of electronic shells OEOJ, as well as of Al clusters, has been obtained [H]. 
In addition, it has been proved that supershells occur naturally in this model [3[7j, which 
is characterized by only one parameter (the deformation parameter q = e T , with r being 
real) in addition to the overall scale. 

An alternative approach to shell structure coming from mathematical physics has been 
introduced long before the discovery of metal clusters [21, by Balian and Bloch [22l, in 
the framework of the theory of periodic orbits and classical quantization conditions |23j. 
In this approach the valence electrons of the metal cluster are supposed to move in a 
spherical cavity (representing the effects of the mean field). Supershells then occur from 
the superposition of closed classical stationary orbits [2H] of the electrons, the superposition 
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of triangular and square orbits being the simplest example The theory of Balian and 
Bloch leads to very specific predictions about the behaviour of magic numbers, which can 

1/3 

serve as a test for alternative approaches. For example, the plot of N i vs. i, where iVj 
are the magic numbers of clusters of a specific alkali metal and % their index [i = 1, 2, 3, 
. . . ), should be a straight line with a slope of 0.61 jU |21j. 

In this paper we confront the results provided by the 3D g-HO to the restrictions imposed 
by the theory of Balian and Bloch. It turns out that agreement between the predictions of 
the two theories can be established by introducing in the Hamiltonian of the 3D g-HO an 
additional term, which is characterized by the same symmetry as the original Hamiltonian. 
Furthermore, it can be seen that this additional term occurs naturally through a variational 
procedure, analogous to the one used for obtaining the Morse oscillator from the usual 
harmonic oscillator. 

In Section 2 of this paper a brief account of the 3D g-HO is given, while in Section 3 
the application of the periodic orbit theory of Balian and Bloch to metal clusters is briefly 
described. The predictions of the 3D g-HO are compared to the restrictions imposed by the 
theory of Balian and Bloch in Section 4, while in Section 5 a modified Hamiltonian for the 
3D g-HO is introduced, allowing for full agreement with the theory of Balian and Bloch. 
Numerical details concerning the modified 3D g-HO Hamiltonian are given in Section 6, 
together with a study of supershells in the framework of this Hamiltonian. In Section 7 
a variational method leading from the usual harmonic oscillator to the Morse oscillator is 
introduced, while in Section 8 this method is applied for deriving the modified Hamiltonian 
introduced in Section 5 from the original 3D g-HO Hamiltonian. Finally in Section 9 a 
discussion of the present results and plans for future work are given. 



2 The 3-dimensional g-deformed harmonic oscillator 
(3D g-HO) 

The space of the 3-dimensional g-deformed harmonic oscillator consists of the completely 
symmetric irreducible representations of the quantum algebra u 9 (3) O 113 El- In this 
space a deformed angular momentum algebra, so 9 (3), can be defined [H3]- The Hamiltonian 
of the 3D g-HO is defined so that it satisfies the following requirements: 

a) It is an so g (3) scalar, i.e. the energy is simultaneously measurable with the g-deformed 
angular momentum related to the algebra so g (3) and its z-projection. 

b) It conserves the number of bosons, in terms of which the quantum algebras u 9 (3) and 
so 9 (3) are realized. 

c) In the limit q — > 1 it is in agreement with the Hamiltonian of the usual 3-dimensional 
harmonic oscillator. 

It has been proved ^Uj that a Hamiltonian of the 3D g-HO satisfying the above require- 
ments takes the form 



H q = hu [N}q 



N+i Q(Q - Q~ 



[2] 




where N is the number operator and Cf® is the second order Casimir operator of the algebra 
so g (3), while 

is the definition of g-numbers and g-operators. In the framework of quantum algebras in 
general, the deformation parameter q can either be real (q = e T , with r being real) or a 
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phase factor (g = e lT , with r being real). In the present case, however, the parameter q is 
restricted to obtain real values only, since this is required in the derivation leading to the 
Hamiltonian of Eq. (1). (See the comment following Eq. (9) in Ref. ^Bj for more details.) 
The energy eigenvalues of the 3D g-HO are given by ^U] 

E q (n,l) = Hu; l[n}q n+1 - 9(g ~ 2] 9 + 1] j = hu e g (n,l), (3) 

where n is the number of vibrational quanta and I is the eigenvalue of the angular momen- 
tum, obtaining the values I — n, n — 2, . . . , or 1, and 

e q (n,l) = [n]g" +1 - g(g [I] [I + 1]. (4) 

In the limit of q — > 1 one obtains \im q ^iE q (n, I) = huofi, which coincides with the 
classical result. 

For small values of the deformation parameter r (where q = e T ) one can expand Eq. (3) 
in powers of r obtaining [10 i 

E q (n, I) = hu n — Tiloqt {1(1 + 1) — n(n + 1)} 

- hu T 2 + 1) - ^n(n + l)(2n + 1)J + C(r 3 ). (5) 

The last expression to leading order bears great similarity to the modified harmonic 
oscillator suggested by Nilsson |2*d1 l2*7j (with the spin-orbit term omitted) . Comparisons 
between the predictions of the 3D g-HO and Nilsson's modified oscillator for the magic 
numbers and supershells of metal clusters have been given in El E|- One of the main 
differences between the two models is that the term rn(n + l) in the eigenvalues of the 
3D g-HO given in Eq. (5) is occuring as a consequence of the overall symmetry, while the 
corresponding term /x'n(n + 3)/2 in the eigenvalues of Nilsson's modified oscillator (see, for 
example, Eq. (3) of Ref. |7j) is put in "by hand". It should also be noticed that the use of 
g-deformations does not mean that an additional parameter is introduced in the theory. In 
both the 3D g-HO and Nilsson's modified oscillator only one parameter appears (in addition 
to the overall scale). In the 3D g-HO case this is the deformation parameter r, while in 
Nilsson's modified oscillator it is the parameter // (see, for example, Eq. (3) of Ref. [?]). 



3 The theory of Balian and Bloch 

The theory of Balian and Bloch j22j is an example of periodic orbit theory [22] , developed as 
a semiclassical bridge between quantum mechanics and classical mechanics. In the approach 
of Balian and Bloch, shell effects in metal clusters are studied by considering the valence 
electrons moving on straight lines within a smooth sphere and being reflected on the inner 
surface of the sphere. Periodic orbits in this case correspond to various polygons with three 
or more corners, the triangular and the square orbits being the simplest examples jipj. 

In the theory of Balian and Bloch a cluster of iV alkali atoms can be studied [2Tj by 
considering the N valence electrons moving in a sphere with radius 

R = r s N 1/3 , (6) 
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where r$ is the Wigner-Seitz radius. The lengths of the triangular and the square orbits 
are then respectively 

L 3 = 3V3TSN 1 / 3 , (7) 

L 4 = 4V2r s N 1/3 . (8) 

Assuming that electrons move with the Fermi velocity vp, the beating pattern created by 
the superposition of the two orbits leads to supershell structure. It turns out that the magic 
numbers (shell closures) iVj are related to the index % counting their number (i = 1, 2, 3, 
• • • ) by m 

h 2 

Nl /3 = — -= = i = 0.605i, (9) 

mv F r s 3^ + 4^ 

where m is the electron mass and the fact that Vprg is constant for all metals has been 
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used. Therefore the plot of N i ' vs. i should be a straight line with a slope of 0.61 . 
Furthermore, a phase shift of a half unit of the index i should occur in the node region of 
the beat pattern, i.e., when passing from a supershell to the next one jUI21]. In addition, 
within each supershell the shell closures should appear at equidistant positions (i.e., they 

should exhibit a periodicity) when plotted vs. N 1 ^ 3 j3] , where N is the number of particles 
(valence electrons in the present case). 

Comparisons to experimental sets of magic numbers have shown 0] that a slope of 0.61 
is obtained in the cases of Na and Li clusters, revealing that the triangular and square 
orbits represent a good approximation in these cases. For Al clusters, though, a slope of 
0.32 has been obtained |4|, indicating that in this case, if the theory is applicable, more 
complicated orbits enter [28] . 



4 Comparing the predictions of the 3D g-HO and of 
the Balian— Bloch theory 

As we have seen in the previous section, a hallmark of the theory of Balian and Bloch is 
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that the plot of N { vs. i, where Ni are the magic numbers and i the index counting them 
(i — 1, 2, 3, . . . ), should be a straight line having a slope of 0.61 in the case of alkali metals, 
while in the case of Al clusters the slope should be 0.32 0j. 

In order to compare the predictions of the 3D g-HO to the restrictions imposed by the 
theory of Balian and Bloch, we plot in Fig. 1(a) (line labelled by e = 0.0) the magic 
numbers obtained from the 3D g-HO for r = 0.038, the parameter value found appropriate 
in Ref. |5| for reproducing the magic numbers of alkali clusters (up to 1500 atoms, which 
is the limit of validity of theories based on the filling of electronic shells [T9"l l20|). Numbers 
considered as magic, listed in Table 1, correspond to energy gaps larger than 5 = 0.38, as 
in Ref. [Sj, with hujQ = 1. In the same figure, a straight line with a slope of 0.61 appears. 
It is clear that the 3D g-HO magic numbers follow the straight line up to % = 14 quite well, 
while beyond this point the predictions of the 3D g-HO are clearly lower than the straight 
line, indicating that "too many" magic numbers are produced by the model in this region. 

The same conclusion is arrived at by looking at Fig. 1(b), where the predictions of the 
3D g-HO for r = 0.050, the parameter value found in Ref. [Sj appropriate for reproducing 
the magic numbers of Al clusters, are reported (line labelled by e = 0.0). Again magic 
numbers, listed in Table 1, correspond to energy gaps larger than 5 = 0.38 (with fkvo = 1). 
Small magic numbers (below 186, i.e., below % = 9) are not shown, since it is known that 
small magic numbers in Al clusters cannot be explained by models based on the filling of 
electronic shells, because of the symmetry breaking caused by the ionic lattice [213 , while 
for large magic numbers this problem does not exist. In the plot a straight line with a slope 
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of 0.32, which is expected to be appropriate for Al clusters jl], as mentioned in the previous 
section, is also seen. It is clear that the predictions of the 3D g-HO follow a line parallel to 
the one with slope 0.32 roughly up to i = 23, while beyond this point the slope is gradually 
reduced, indicating that the model predicts "too many" magic numbers. 



5 A modified Hamiltonian for the 3D g-HO 

The discrepancy between the 3D g-HO and the theory of Balian and Bloch can be lifted by 
considering the Hamiltonian 

H' q = H q -eHl (10) 

with eigenvalues 

E' q (n,l)=E q (n,l)-eE 2 q (n,l), (11) 

where e is a small real positive constant. Justification for this choice will be given in Sections 
7 and 8 through a variational procedure. For the moment the following comments suffice: 

a) It is clear that H' is a function of H q , which is by construction an so 9 (3) scalar, as 

mentioned in Section 2. Therefore H' is also an so g (3) scalar. 

b) The energy eigenvalues of the new Hamiltonian can be written in the form 

E' (n, I) = Hujq (1 — ehuJoe q (n, /)) e q (n, I) = hu(n, l)e q (n, I), (12) 

where 

L)(n, I) — ojq (1 — efiujQe q {n, I)) (13) 

is a variable frequency, depending on the quantum numbers n, I, and on the small parameter 
e. We shall call this oscillator the Variable Frequency Oscillator (VFO) corresponding to 
the 3D g-HO, a term for which justification will be provided in Sections 7 and 8. 

The magic numbers provided by the VFO for a few appropriate values of e in the case 
of r = 0.038, which is relevant for alkali clusters (5., are shown in Fig. 1(a) and listed in 
Table 1. Details of the calculation will be given in Section 6. Once more magic numbers 
are separated by gaps larger than 5 = 0.38, while Tiujq = 1. It is clear that the predictions 
of the VFO roughly follow the straight line with a slope of 0.61 even for large values of i, 
thus overcoming the difficulties faced by the 3D g-HO. 

A similar picture is obtained for r = 0.050 (and 5 = 0.38, with huoo = 1), which is 
appropriate for Al clusters 0. For a few appropriate values of e, shown in Fig. 1(b) and 
listed in Table 1, the predictions of the VFO roughly follow a straight line with a slope of 
0.32, as they should, according to the previous section. 

We therefore conclude that the addition of the second term in the Hamiltonian of the 
3D g-HO, leading to the VFO, makes the predictions of the 3D g-HO compatible with the 
predictions of the theory of Balian and Bloch. 



6 Numerical details 

In this section the calculations leading to the results reported in Section 5 will be described. 
Throughout this paper we put %ujq = 1 for simplicity. 

An important difference between the 3D g-HO and the VFO of Eq. (fTU|) lies in the way 
truncations of the spectrum are made. The following comments apply: 

a) In the 3D g-HO the level with I = n always lies lowest in energy within each shell, 
the level with I = n — 2 lies immediately above it, and so on. Therefore stopping the level 
scheme at the / = n level of a given shell and taking into account all levels with lower n 
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(i.e., all levels of the shells lying below the given one), one makes sure that all levels up to 
the given level have been included [Zj. 

b) In the case of the VFO of Eq. (JTUJl the following picture occurs: For a given (small) 
value of e the first several shells exhibit the same behaviour as in the case of the 3D (/-HO, 
i.e., the level with I = n lies lowest in energy within each shell, the level with / = n — 2 lies 
immediately above it, and so on. As the shell number is increasing, however, an inversion 
of the order of the levels occurs, with the levels with I = (for n being even) or I = 1 
(for n being odd) lying lowest in energy within the shell. It is then clear that beyond this 
inversion point truncation at a given shell should be made at the level with I = (for n 
being even) or at I — 1 (for n being odd). 

The magic numbers obtained for the values of r and e used in the previous section are 
given in Table 1. In each case the maximum value of n, n max , included in the calculation 
is indicated. Care has been taken that no inversion of the order of the levels, of the type 
described in comment b) above, occurs for the values of n included in the calculation. 
Therefore in all cases truncation is made at the level with n = n max and I = n max . The 
total number of levels up to the truncation point, N max , is also shown in Table 1. The 
following remarks are now in place: 

a) In the cases considered here and up to the truncation point, for a given value of r 
the order of the levels is not modified as e is changed. The only modification occuring is 
that the spectrum gets "squeezed" as e increases. 

b) As a result of a), the magic numbers reported in Table 1 present the following feature. 
For each value of r and for e = 0.0 the largest number of magic numbers appears. As e 
increases, some of the magic numbers cease to be magic any more, since the "squeezing" of 
the spectrum brings the levels closer to each other. For each value of e the magic numbers 
occuring are a subset of the magic numbers occuring for lower values of e (with the same 
r). No new magic numbers appear, within the limits considered here, as e increases. 

It is interesting to examine at this point what the influence of the additional term to the 
appearance of supershells is. For this purpose we are going to use the procedure employed 
by Nishioka et al. [30, 31J. For a given number of particles N the single particle energies 
Ej(n, I) of the N occupied states are summed up 

JV 

E(N) = Y,E 3 (n,l). (14) 

3=1 

This sum is then divided into two parts: A smooth average part E av and a shell part E s h e u, 
which will exhibit the supeshell structure 

E(N) = E av (N)+E shell (N). (15) 

For the average part of the total energy a Liquid Drop Model expansion is used [Zj 

E av (N) = ai N 1/3 + a 2 N 2/3 + a 3 N + a 4 iV 4/3 + a 5 N 5/3 + a 6 N 2 . (16) 

The parameters of the fits occuring in the cases considered here are shown in Table 2, 
together with the number of levels, N max , included in the fit and the rms deviation a. We 
remark that for a given value of r the parameters change smoothly for the nonzero values 
of e, while the case with e = 0.0 is characterized by quite different values of the parameters, 
but also by a higher rms deviation a. The addition of the second term in Eq. (jlOj) improves 
the agreement of the average part of the total energy to the Liquid Drop Model expansion, 
thus resulting in lower rms deviations a. 

The procedure of the calculation was as follows: First the summations described by Eq. 
()14|) . resulting in the total energy E(N) for each particle number N, have been performed. 
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Subsequently, in order to reduce the size of the calculation approximately by a factor of 
10, the average E(N) was calculated every 11 points (i.e., for N = 6, 17, 28, ... ) up to 
the cutoff point which is reported in Table 2 as N max . These averaged values of E(N) 
were subsequently fitted by the expansion of Eq. (fTBj). resulting in the determination of 
E av (N) at these points. Finally E sheU (N) has been obtained at these points as the difference 
E(N) - E av (N) and plotted in Figs. 2 and 3. 

The shell energy, Egheii, is plotted vs. the particle number, N, for r = 0.038, which is 
appropriate for alkali clusters [2], in Fig. 2. In Fig. 2(a) the predictions of the original 3D 
g-HO (with e = 0) are shown, while in Fig. 2(b) the results of the VFO with e = 0.006 are 
depicted. The VFO with e = 0.007, 0.008 gives results which look almost identical with 
Fig. 2(b) and therefore are not shown for brevity. The similarities between Figs. 2(a) and 
2(b) are clear. Not only the supershell appears in both cases around N = 1000, as it is 
expected for Na clusters [HOI EH], but in addition even the maxima and minima of the shell 
energy appear at the same particle numbers and have roughly the same magnitude. Even 
the local maxima and minima present striking similarities. These results corroborate the 
remarks made above, namely that the addition of the second term to the Hamiltonian of 
Eq. (fTUj) does not influence the order of the energy levels, the main effect of the second 
term being the gradual "squeezing" of the spectrum as energy increases. Of course this 
conclusion is valid only within the region of particle numbers studied and for small values 
of e, like the ones used here. 

Similar results are obtained in Fig. 3 for the case of r = 0.050, which is appropriate 
for Al clusters 5 . Besides the e = case, shown in Fig. 3(a), the results corresponding to 
e = 0.0050 are shown in Fig. 3(b), since the cases with e = 0.0053, 0.0055 provide results 
almost identical with the ones shown in Fig. 3(b). In all cases there is some evidence 
for a supershell below N = 1000, although its appearance is not as clear as in the case 
of Fig. 2. The appearance of a supershell in this region is in agreement with the results 
of more sophisticated calculations, as, for example, spherical jellium model predictions in 
Local Density Approximation [32J, but it is not in good agreement with experiment, where 
no evidence for supershell in Al clusters exists in this region [Hill EH]- An advantage of the 
VFO in comparison to earlier calculations [H2] is that at least it can reproduce the slope of 
0.32 in Fig. 1(b), something which is not occuring in spherical jellium model calculations, 
although it occurs experimentally jUE3 Ell- 
in the theory of Balian and Bloch, as mentioned above, a phase shift by a half unit of 
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the running index i should be observed in the plot of the magic numbers N i vs. i when 
passing from a supershell to the next [UI23]. 111 Fig- l( a ) such a shift is seen quite clearly 
around N = 1000 in the cases of e = 0.006, 0.007, 0.008, while no clear shift of this type is 
seen in Fig. 1(b). 

Furthermore, the theory of Balian and Bloch, as mentioned above, predicts that within 
each supershell the minima of the shell energy, E sheU) should appear at equidistant positions 
(i.e. they should exhibit a periodicity) when plotted vs. V 1 / 3 . The change in the periodicity 
when passing from the first supershell to the second one is clear in Figs. 2 and 3 (although 

in these cases, for reasons of clarity, E s h e u is plotted vs. N and not vs. iV 1//3 ). The fact 
that the predictions of the VFO corresponding to the 3D g-HO approximately exhibit the 
right periodicity features does not come as a surprise, since the 3D g-HO is known to show 
this feature [HI El, while, as we have seen above, the addition of the second term in Eq. (fTUjl 
does not influence the position of the minima. 

In conclusion, the VFO corresponding to the 3D g-HO is able to reproduce the right 
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slope in the N t vs. i plot in both the alkali and Al clusters. In addition it predicts correctly 
the first supershell in alkali clusters, while in Al clusters its prediction for a supershell is in 
rough agreement with results of spherical jellium models but not with experiment. 
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7 A variational method 



In nuclear physics it is well known that nuclear spectra can be described very accurately in 
terms of the Variable Moment of Inertia (VMI) model [HE]- in this model the energy levels 
are given by 

E(J) = + \C{Q{J) - 6(0)} 2 , (17) 

where J is the angular momentum and Q(J) is the moment of inertia, which is supposed 
to be a function of the angular momentum. C and 0(0) are free parameters, the latter 
representing the ground state moment of inertia. It is clear that the VMI formula is a 
generalization of the rigid rotator formula 
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in which the moment of inertia is assumed to be constant. The rigid rotator formula is 
known to fail beyond the first few levels of a rotational nucleus, since the experimental 
levels appear "squeezed" in comparison to the rigid rotator predictions. This difficulty is 
overcome in the framework of the VMI model by determining the moment of inertia for 
each value of the angular momentum J through a minimization of the energy with respect 
to the moment of inertia for given angular momentum 

aeM |j = ' (19) 

This variational condition leads to a cubic equation for O(J), which turns out to have only 
one real solution [33], corresponding to the appropriate value of the moment of inertia for 
the given value of the angular momentum. The second term in Eq. (jl7J) is justified by the 
well known fact that many perturbing potentials near their origin can be approximated by 
a harmonic oscillator potential. 

Following the same reasoning, it is interesting to examine what happens to the usual 
harmonic oscillator if, by analogy, one allows the angular frequency to be a function of the 
quantum number n. The energy will then read 

E(n) = hu(n) (n + + -C{u(n) - cu(0)} 2 , (20) 

where C and u;(0) are free parameters, the latter corresponding to the ground state angular 
frequency. The variational principle in this case should correspond to the minimization of 
the energy with respect to the angular frequency for constant value of the quantum number 

n 

dE(n) 



duj{n) 

It is clear that this condition leads to 
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h (n + + C{u{n) - lj(0)} = u(n) = u(0) - ^ (n + ^\ . (22) 
Substituting this result in Eq. (|2T)j) we obtain 

E(n) = MO) [n + £) - ~ [n + , (23) 
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which is reminiscent of the spectrum of the Morse potential . 

Indeed, solving the Schrodinger equation for the Morse potential [0E1 EH IHHj 



V(x) = D{1 - e- ax )\ (24) 

one obtains the energy spectrum 




2" 



(25) 



where 

I ha , 
x e = — ; , 26 



and 




(27) 



We therefore conclude that by allowing the angular frequency of the simple harmonic 
oscillator to vary with the quantum number n, we obtain the spectrum of the Morse oscil- 
lator. Again the second term in Eq. (J20|) is in agreement to the fact that most perturbing 
potentials near their origin can be app roximated by the harmonic oscillator potential. We 
shall refer to the oscillator of Eq. (}20j) as the Variable Frequency Oscillator (VFO). 



8 Derivation of the modified Hamiltonian for the 3D 
g-HO 

through a variational method 

The idea leading to the VFO of the previous section can be appropriately generalized in 
the case of the 3D q-HO. In this case we consider the energy expression 

E' q (n,l) = huj(n,l) {[n]^ 1 - q{q ~ 2 ^ ] [l][l + 1]} + \c{u{n,l) -c(0,0)} 2 

= huj(n, l)e g (n, I) + -C{u(n, I) - u(0, 0)} 2 , (28) 

where the angular frequency u(n,l) depends on the quantum numbers n and I, while C 
and u(0, 0) are real positive constants, the latter corresponding to the ground state angular 
frequency, since the ground state of the 3D g-HO is characterized by n = and / = 0. It 
is thus clear that u(0,0) corresponds to ojq appearing in Eq. (jSJ), i.e., cu(0,0) = u . The 
variational condition in the present case should read 

idJ k ' = - (29) 

In other words, the energy is minimized with respect to the angular frequency for constant 
values of the quantum numbers n and I. The variational condition leads to the equation 

he q (n, I) + C{lu(u, I) - u(0, 0)} = =>• w(n, I) = u(0, 0) - ^e g (n, I) 
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"(0, 0) - M [n]<f +1 - q{q [2 p [I] [/ + !]}. (30) 



in 2 

E' q (n, I) = MO, 0)e,(n, Z) - - -ej(n, /), (31) 



Substituting this result in Eq. (jzHj) one then obtains 



which is the same as Eq. (fTT|) . with 



since u = c<j(0, 0), as mentioned above. 

9 Discussion 

In this paper we have attempted a comparison of the predictions for the shell structure 
of metal clusters of the 3D g-HO model to the ones of the periodic orbit theory of Balian 
and Bloch. It turns out that the predictions of the 3D g-HO for the magic numbers of 
metal clusters can be made compatible with the predictions of the theory of Balian and 
Bloch by adding to the 3D g-HO a symmetry-preserving correction term reminiscent of the 
anharmonicity term in the spectrum of the Morse potential, while this addition does not 
influence the predictions for the supershells. This extended expression for the 3D g-HO can 
be justified through a variational method, similar to the one used in the Variable Moment 
of Inertia (VMI) model of nuclear physics, leading to the concept of the Variable Frequency 
Oscillator (VFO), which gives promise of wider applicability. 

Acknowledgements 

One of the authors (PPR) acknowledges support from the Bulgarian Ministry of Science 
and Education under contracts $-415 and $-547. 



References 

[1] F. Iachello and A. Arima, The Interacting Boson Model (Cambridge University Press, 
Cambridge, 1987). 

[2] F. Iachello and R. D. Levine, Algebraic Theory of Molecules (Oxford University Press, 
Oxford, 1995). 

[3] W. A. de Heer, Rev. Mod. Phys. 65, 611-676 (1993). 

[4] M. Brack, Rev. Mod. Phys. 65, 677-732 (1993). 

[5] D. Bonatsos, N. Karoussos, D. Lenis, P. P. Raychev, R. P. Roussev and P. A. Terziev, 
Phys. Rev. A 62, 013203 [13 pages] (2000). 

[6] D. Bonatsos, D. Lenis, P. P. Raychev and P. A. Terziev, Int. J. Quant. Chem. 89, 
299-312 (2002). 



10 



[7] D. Bonatsos, D. Lenis, P. P. Raychev and P. A. Terziev, Phys. Rev. A 65, 033203 [12 
pages] (2002). 

[8] A. I. Kuleff, J. Maruani and P. P. Raychev, Adv. Quant. Chem. 40, 279-304 (2001). 

[9] D. Bonatsos, A. I. Kuleff, J. Maruani, P. P. Raychev and P. A. Terziev, Int. J. Quant. 
Chem. 89, 377-388 (2002). 

[10] P. P. Raychev, R. P. Roussev, N. Lo Iudice and P. A. Terziev, J. Phys. G: Nucl. Part. 
Phys. 24, 1931-1943 (1998). 

[11] D. Bonatsos and C. Daskaloyannis, Prog. Part. Nucl. Phys. 43, 537-618 (1999). 

[12] Yu. F. Smirnov, V. N. Tolstoy and Yu. I. Kharitonov, Yad. Fiz. 54, 721-736 (1991) 
[Sov. J. Nucl. Phys. 54, 437-445 (1991)]. 

[13] J. Van der Jeugt, J. Phys. A: Math. Gen. 25, L213-L218 (1992). 

[14] J. Van der Jeugt, J. Math. Phys. 34, 1799-1806 (1993). 

[15] C. Quesne, Phys. Lett. 5 304, 81-88 (1993). 

[16] P. P. Raychev, R. P. Roussev, P. A. Terziev, D. Bonatsos and N. Lo Iudice, J. Phys. 
A: Math. Gen. 29, 6939-6949 (1996). 

[17] Yu. F. Smirnov, V. N. Tolstoy and Yu. I. Kharitonov, Yad. Fiz. 53, 959-980 (1991) 
[Sov. J. Nucl. Phys. 53, 593-605 (1991)]. 

[18] Yu. F. Smirnov, V. N. Tolstoy and Yu. I. Kharitonov, Yad. Fiz. 56, 223-244 (1993) 
[Phys. At. Nucl. 56, 690-700 (1993)]. 

[19] T. P. Martin, T. Bergmann, H. Gohlich and T. Lange, Chem. Phys. Lett. 172, 209-213 
(1990). 

[20] K. Clemenger, Phys. Rev. B 44, 12991-13001 (1991). 

[21] W. D. Knight, K. Clemenger, W. A. de Heer, W. A. Saunders, M. Y. Chou and M. L. 
Cohen, Phys. Rev. Lett. 52, 2141-2143 (1984). 

[22] R. Balian and C. Bloch, Ann. Phys. (N.Y.) 69, 76-160 (1972). 

[23] M. C. Gutzwiller, J. Math. Phys. 12, 343-358 (1971). 

[24] C. Brechignac, Ph. Cahuzac, F. Carlier, M. de Frutos and J. Ph. Roux, Phys. Rev. B 
47, 2271-2277 (1993). 

[25] P. Morse, Phys. Rev. 34, 57-64 (1929). 

[26] S. G. Nilsson, Mat. Fys. Medd. K. Dan. Vidensk. Selsk. 29, 16 (1955). 

[27] S. G. Nilsson and I. Ragnarsson, Shapes and Shells in Nuclear Structure (Cambridge 
University Press, Cambridge, 1995). 

[28] J. Lerme, M. Pellarin, J. L. Vialle, B. Baguenard and M. Broyer, Phys. Rev. Lett. 68, 
2818-2821 (1992). 



11 



[30 
[31 
[32 

[33 

[34; 

[35 
[36 

[37; 

[38 



J. L. Persson, R. L. Whetten, H. P. Cheng and R. S. Berry, Chem. Phys. Lett. 186, 
215-222 (1991). 

H. Nishioka, K. Hansen and B. R. Mottelson, Phys. Rev. B 42, 9377-9386 (1990). 

H. Nishioka, Z. Phys. D 19, 19-23 (1991). 

0. Genzken, M. Brack, E. Chabanat and J. Meyer, Ber. Bunsenges. Phys. Chem. 96, 
1217-1220 (1992). 

M. Pellarin, J. Lerme, B. Baguenard, M. Broyer and J. L. Vialle, Ber. Bunsenges. 
Phys. Chem. 96, 1212-1215 (1992). 

M. Pellarin, B. Baguenard, M. Broyer, J. Lerme and J. L. Vialle, J. Chem. Phys. 98, 
944-950 (1993). 

M. A. J. Mariscotti, G. Scharff-Goldhaber and B. Buck, Phys. Rev. 178, 1864-1887 
(1969). 

S. Fliigge, Practical Quantum Mechanics (Springer, Berlin, 1974). 

I. L. Cooper, Chem. Phys. 112, 67-75 (1987). 

1. L. Cooper, J. Phys. A: Math. Gen. 25, 1671-1683 (1992). 



Figure captions 

Fig. 1 Cubic roots of the magic numbers iVj plotted vs. the running index i counting 
them. The magic numbers are the ones listed in Table 1. (a) r = 0.038, (b) r = 0.050. 

Fig. 2 Shell part (E s h e u) of the total energy [in units of huo, see Eqs. (JHJ) and fTTf ] 
vs. the number of particles N, in the case of r = 0.038, (a) for the 3D q-B.0, (b) for a 
corresponding VFO. The values of the dimensionless parameters r and e are listed in Table 
2, together with the details of the calculation. See Section 6 for further discussion. 

Fig. 3 Same as Fig. 2, but for r = 0.050. 
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Table 1: Magic numbers (corresponding to gaps larger than 5 = 0.38, with hooo = 1) 
produced by the 3D g-HO [Eq. (j3J)] (cases with e = 0.0) and the corresponding VFO [Eq. 
([lip] for different values of the parameters r and e. n max is the maximum value of the 
quantum number n included in the calculation, while the highest level taken into account 
is the one with n = n max and I = n max , corresponding to the reported total number of 
particles N max . i is a running index counting the magic numbers. See Section 6 for further 
discussion, 



T 


0.038 


0.038 


0.038 


0.038 


0.050 


0.050 


0.050 


0.050 


e 


0.0 


0.006 


0.007 


0.008 


0.0 


0.0050 


0.0053 


0.0055 


^max 


26 


26 


26 


25 


26 


26 


25 


25 


N 

1 v max 


4658 


4658 


4658 


4154 


4778 


4778 


4258 


4258 


i 
1 


2 


2 


2 


2 


2 


2 


2 


2 


2 


8 


8 


8 


8 


8 


8 


8 


8 


3 


20 


20 


20 


20 


20 


20 


20 


20 


4 


34 


34 


40 


40 


34 


34 


34 


34 


5 


40 


40 


58 


58 


40 


40 


40 


40 


6 


58 


58 


92 


92 


58 


58 


58 


58 


7 


92 


92 


138 


138 


92 


92 


92 


92 


8 


138 


138 


198 


198 


138 


138 


138 


138 


9 


198 


198 


254 


254 


186 


186 


186 


186 


10 


254 


254 


338 


338 


254 


254 


254 


254 


11 


268 


338 


440 


440 


338 


338 


338 


338 


12 


338 


440 


676 


676 


398 


398 


398 


398 


13 


440 


676 


832 


832 


440 


440 


440 


440 


14 


556 


832 


912 


912 


486 


542 


542 


542 


15 


562 


912 


1012 


1012 


542 


612 


612 


612 


16 


676 


1012 


1100 


1100 


612 


676 


676 


676 


17 


694 


1100 


1206 


1206 


676 


748 


748 


748 


18 


832 


1206 


1660 


1660 


748 


832 


832 


832 


19 


912 


1660 


1760 


1760 


832 


912 


912 


912 


20 


1012 


1760 


2048 


2048 


890 


1006 


1006 


1006 


21 


1100 


2048 


2368 


2368 


912 


1074 


1074 


1074 


22 


1206 


2368 


3028 


3028 


1006 


1100 


1100 


1100 


23 


1284 


3028 


3438 


3438 


1074 


1284 


1284 


1284 


24 


1314 


3438 


3886 


3886 


1100 


1314 


1314 


1410 


25 


1410 


3886 


4374 




1206 


1410 


1410 


1502 


26 


1502 


4052 






1284 


1502 


1502 


1760 


27 


1516 


4374 






1314 


1516 


1760 


2018 


28 


1660 








1410 


1760 


2018 


2048 


29 


1760 








1502 


2018 


2048 


2178 


30 


2018 








1516 


2048 


2178 


2334 
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Table 1: (continued) 



T 


0.038 


0.038 


0.038 


0.038 


0.050 


0.050 


0.050 


0.050 


€ 


0.0 


0.006 


0.007 


0.008 


0.0 


0.0050 


0.0053 


0.0055 




26 


26 


26 


25 


26 


26 


25 


25 




4658 


4658 


4658 


4154 


4778 


4778 


4258 


4258 


i 

31 


2048 








1614 


2178 


2334 


2368 


32 


2178 








1660 


2334 


2368 


2510 


33 


2334 








1734 


2368 


2510 


2672 


34 


2368 








1760 


2510 


2672 


2722 


35 


2654 








1778 


2672 


2722 


3028 


36 


2672 








1940 


2722 


3028 


3050 


37 


2722 








2018 


3028 


3050 


3112 


38 


2796 








2048 


3050 


3112 


3438 


39 


3028 








2178 


3112 


3438 


3464 


40 


3050 








2334 


3438 


3464 


3886 


41 


3190 








2368 


3464 


3886 


3916 


42 


3404 








2510 


3886 


3916 




43 


3438 








2672 


3916 






44 


3464 








2684 


3988 






45 


3610 








2722 


4374 






46 


3848 








2876 


4408 






47 


3886 








3028 








48 


4052 








3050 








49 


4312 








3112 








50 


4326 








3190 








51 


4374 








3244 








52 


4552 








3438 








53 










3464 








54 










3528 








55 










3622 








56 










3680 








57 










3886 








58 










3916 








59 










3988 








60 










4088 








61 










4156 








62 










4374 








63 










4408 








64 










4462 








65 










4488 








66 










4578 








67 










4596 
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Table 2: Parameters used for fitting the average part of the total energy [see Eq. (|16j) ] in the 
case of the 3D g-HO and the corresponding VFO for various values of the parameters r and 
e, corresponding to the cases exhibited in Figs. 2 and 3. The parameters are dimensionless, 
since we have assumed hu>o = 1 [see Eqs. © and (JTTJ)] throughout. The number of particles 
N max included in each calculation and the relevant rms deviation a are also shown. See 
Section 6 for further discussion. 



10 3 r 


10 4 e 


a\ 


a 2 


a 3 


10a4 


10 2 a 5 


10 4 a 6 


Nmcix 


a 


38 





-21.035 


18.295 


-7.295 


19.521 


-6.082 


40.857 


3009 


8.904 


38 


60 


24.756 


-20.883 


5.201 


0.493 


7.946 


-8.993 


3009 


5.758 


38 


70 


32.475 


-27.496 


7.306 


-2.704 


10.297 


-17.326 


3009 


5.297 


38 


80 


39.762 


-33.786 


9.329 


-5.806 


12.597 


-25.559 


3009 


4.834 


50 





-24.946 


24.641 


-10.384 


26.117 


-13.000 


82.558 


2008 


7.328 


50 


50 


14.051 


-13.208 


3.323 


2.409 


7.006 


6.484 


2008 


5.286 


50 


53 


16.795 


-15.857 


4.264 


0.817 


8.320 


1.634 


2008 


5.175 


50 


55 


18.589 


-17.556 


4.862 


-0.188 


9.150 


-1.464 


2008 


5.098 
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